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ABSTRACT 

We present 2.5D radiation-hydrodynamics simulations of the accretion-induced collapse (AIC) of 
white dwarfs, starting from 2D rotational equilibrium configurations of a 1.46-M0 and a 1.92-M0 
model. Electron capture leads to the collapse to nuclear densities of these cores within a few tens 
of milliseconds. The shock generated at bounce moves slowly, but steadily, outwards. Within 50- 
100 ms, the stalled shock breaks out of the white dwarf along the poles. The blast is followed by a 
neutrino-driven wind that develops within the white dwarf, in a cone of ^ 40° opening angle about 
the poles, with a mass loss rate of 5-8x 10~^MQyr~^. The ejecta have an entropy on the order of 
20-50 ke/baryon, and an electron fraction distribution that is bimodal. By the end of the simulations, 
at >600ms after bounce, the explosion energy has reached 3-4x10''^ erg and the total ejecta mass 
has reached a few times 0.001 M©. We estimate the asymptotic explosion energies to be lower than 
10^° erg, significantly lower than those inferred for standard core collapse. The AIC of white dwarfs 
thus represents one instance where a neutrino mechanism leads undoubtedly to a successful, albeit 
weak, explosion. 

We document in detail the numerous effects of the fast rotation of the progenitors: The neutron 
stars are aspherical; the 'V^" and i>f. neutrino luminosities are reduced compared to the Uf, neutrino 
luminosity; the deleptonized region has a butterfly shape; the neutrino flux and electron fraction 
depend strongly upon latitude (d la von Zeipel); and a quasi-Keplerian 0.1-0.5-Mq accretion disk is 
formed. 

Subject headings: hydrodynamics - neutrinos - rotation - stars: neutron - stars: supernovae: general 
- stars: white dwarfs 



1. INTRODUCTION 

Stars can follow a few special evolutionary routes to 
form an unstable Chandrasekhar mass core. A main- 
sequence star of more than ~8 M© evolves to form either 
a degenerate 0/Ne/Mg core (Barkat et al. 1974; Nomoto 
1984,1987; Miyaji & Nomoto 1987) or a degenerate Fe 
core (Woosley & Weaver 1995), which, due to photodis- 
integration of heavy nuclei and/or electron capture, col- 
lapses to form a protoneutron star (PNS). If an explosion 
ensues, the event is associated with a Type II supernovae 
(SN). Less massive stars end their lives as white dwarfs. 
White dwarfs located in a binary system may accrete 
from a companion and achieve the Chandrasekhar mass, 
triggering the thermonuclear runaway of the object and 
leading to Type la SN, leaving no remnant behind. 

However, a third class of objects is expected. The- 
oretically, massive white dwarfs with 0/Ne/Mg cores, 
due to their high central density {>10^° gcm"'^), expe- 
rience rapid electron capture that leads to the collapse 
of the core. This is accretion- induced collapse (AIC), 
an alternative path to stellar disruption through ex- 
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plosive burning, currently associated with Type la SN 
(Nomoto & Kondo 1991). It is presently unclear what 
fraction of all white dwarfs will lead to AICs, but of those 
white dwarfs that evolve to form a Chandrasekhar-mass 
0/Ne/Mg core, all will necessarily undergo core collapse. 
One formation channel is the coalescence of two white 
dwarfs (Mochkovitch & Livio 1989), with either C/0 
or 0/Ne/Mg cores, although few such binary systems 
have yet been observed with a cumulative mass above 
the Chandrasekhar mass. There is still uncertainty as 
to whether such binary systems would not undergo ther- 
monuclear runaway rather than collapse. Since the coa- 
lescence of two white dwarfs requires a shrinking of the 
orbit through gravitational radiation, these systems will 
take many gigayears to coalesce. An alternative forma- 
tion mechanism is via single-degenerate systems, through 
a combination of high original white dwarf mass and 
mass and angular-momentum accretion by mass transfer 
from a (non-degenerate) H/He star (Nomoto & Kondo 
1991). Binary star population synthesis codes predict 
the occurence of the AIC of white dwarfs with a galac- 
tic rate of 8xl0~^ yr~-^ to 8xl0~^ yr~^, depending, 
amongst other things, on the treatment of the common- 
envelope phase and mass transfer (Yungelson & Livio 
1998). The set of parameters leading to a Type la rate 
of 10^'^ yr^^ corresponds to an AIC rate of 5xl0~^ yr~^. 
The observed Type la rate of ^3xl0~'^yr~^ (Madau et 
al. 1998; Blanc et al. 2004; Manucci et al. 2005) would 
imply a galactic AIC rate of 1.5x 10~* yr~^. These rates 
are likely functions of galaxy and metallicity (Yungel- 
son & Livio 2000; Belczynski et al. 2005; Greggio 2005; 
Scannapieco & Bildsten 2005). Based on r-process nu- 
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cleosynthetic yields obtained from previous simulations 
of the AIC of white dwarfs, Fryer et al. (1999) inferred 
rates ranging from ~10~^ to ^--^10"^ yr^^. Overall, AICs 
are not expected to occur more than once per 20-50 stan- 
dard Type la events; because they are intrinsically rarer, 
they remain to be identified and observed in Nature. 

Whatever their origin, their fundamental nature is to 
accrete both mass and angular momentum from a com- 
panion object. Rotation is, therefore, a key physical 
component. Prior to core collapse, differential rotation 
acts as a stabilizing agent for shell burning by widening 
its spatial extent through enhanced mixing and reducing 
the envelope density through centrifugal support (Yoon 
& Langer 2004) . Mass accretion also leads to an increase 
of the central density, which may rise up to a few 10^*^ 
gcm~'^, establishing suitable conditions for efhcient elec- 
tron capture on Mg/Ne nuclei. By the time of core col- 
lapse and depending on the evolutionary path followed, 
such white dwarfs may cover a range of masses from 
>1.35M0up to ~2M0 (2.7Mq) in the case of a non- 
degenerate (degenerate) companion, potentially well in 
excess of the standard Chandrasekhar limit, and possess- 
ing an initial rotational energy up to ^--^10% of their gravi- 
tational binding energies. At such values, the centrifugal 
potential leads to a deformation of the white dwarf from 
spherical symmetry, equipotentials and isopressure sur- 
faces adopting a peanut-like shape in cross section for 
the largest rotation rates. Such structures are obtained 
in the 2D differentially rotating equilibrium white dwarf 
models constructed by Yoon & Langer (2005, YL05; see 
also Liu & Lindblom 2001). 

In the past, the collapse of 0/Ne/Mg cores originat- 
ing from stars in the 8-10 M0 range has been studied in 
ID by Baron et al. (1987ab), Mayle & Wilson (1988), 
and Woosley & Baron (1992), who showed that the shock 
generated at core bounce stalls rather than leading to a 
prompt explosion. However, Hillebrandt et al. (1984) 
and Mayle & Wilson (1988) obtained delayed explosions, 
the former after 20-30 ms, and the latter after ~200ms, 
supposedly driven by neutrino energy deposition behind 
the stalled shock. Woosley & Baron (1992) found the 
emergence of a sustained neutrino-driven wind, with a 
mass loss rate of 0.005 M© s~^ and with an ejecta elec- 
tron fraction of the order of 0.45, showing promise, mod- 
ulo uncertainties, for a contribution to the enrichment of 
the ISM in r-process elements. Using 1D/2D SPH sim- 
ulations. Fryer et al. (1999) reproduced the simulations 
of Woosley & Baron (1992), confirmed that the shock 
stalls due to the copious neutrino losses associated with 
core bounce, did not find prompt explosion, and, depend- 
ing on the equation of state (EOS) employed, observed 
a delayed explosion. They focused mostly on the early 
phase, prior to the neutrino-driven wind, and found an 
ejected mass of low-l^ material of ~0.05 Mq , depend- 
ing on adopted model assumptions. Their 2D simula- 
tions with solid-body rotation showed similar properties 
to their ID equivalents, the authors attributing the small 
differences to the different grid resolution. This may 
partly stem from the essentially spherical explosion trig- 
gered just <100ms after core bounce, ejecting the fast- 
rotating material in the outer mantle. Consequently, at 
the end of their 2D simulations, they obtain very slow ro- 
tation rates for the PNS, i.e., of ~ls. Recently, Kitaura 
et al. (2005) re-inspected the collapse of the progenitor 



model used by Hillebrandt et al. (1984). They confirmed 
again the by-now well-accepted idea that no prompt ex- 
plosion occurs, but instead obtain a successful, though 
sub-energetic, delayed explosion in spherical symmetry, 
powered by neutrino heating and a neutrino-driven wind 
that sets in '^200ms after bounce. 

A primary motivation for this work is to improve 
upon these former investigations that assumed one- 
dimensionality, sphericity, and/or zero- rotation, and 
start instead from the more physically-consistent 2D 
models of YL05, thereby fully accounting for the effects 
of rotation on the collapse, bounce, and post-bounce evo- 
lution of the white dwarf core and envelope, as well as 
for the strong asphericity of the progenitor. Our study 
uses VULCAN/2D (Livne et al. 2004; Walder et al. 
2005) to perform 2D Multi-Group Flux Limited Dif- 
fusion (MGFLD) radiation hydrodynamics simulations. 
As we will demonstrate, rotation plays a major role in 
the post-bounce evolution, making ID investigations of 
such objects of limited utility. By carrying out the sim- 
ulations from ~30ms before bounce to >600ms after 
bounce, we capture a wide range of physical processes, 
including the establishment of a strong, fast, and aspher- 
ical neutrino-driven wind. We model the centrifugally- 
supported equatorial regions and the large angular mo- 
mentum budget leading to the formation of a sizable 
accretion disk. Moreover, an Eulerian investigation is 
better-suited than a Lagrangean approach to explore the 
neutrino-driven wind that develops after >200 ms. Fi- 
nally, despite the relative scarcity of AIC in Nature, these 
simulations represent interesting examples for the forma- 
tion of disks around neutron stars. 

The main findings of this work are the following: We 
find that the AIC of white dwarfs forms ~1.4-M0 neutron 
stars, expelling a modest mass of a few 10"'^ Mq mostly 
through a neutrino-driven wind that develops >200ms 
after bounce, and that they lead to very modest explo- 
sion energies of 5-lOx 10''^erg. Accounting for the rota- 
tion and the asphericity of the progenitor white dwarfs 
reveals a wealth of phenomena. The shock wave gen- 
erated at core bounce emerges through the poles rather 
than the equator, and it is in this excavated polar region 
that the neutrino-driven wind develops. The strong as- 
phericity of the newly- formed protoneutron stars leads to 
a latitudinally-dependent neutrino flux, while the effects 
of rotation modify the relative flux magnitude of differ- 
ent neutrino flavors. Besides mass accretion by both the 
neutron star and mass ejection by the initial blast and 
the subsequent wind, we find a sizable component that 
survives and resides in a quasi-Keplerian disk, which ob- 
structs the wind fiow at low latitudes. This disk will 
be accreted by the neutron star only on longer, viscous 
timescales. 

In the next two sections, we present the two se- 
lected progenitor models in more detail; we also discuss 
the radiation- hydrodynamics code VULCAN/2D and the 
various assumptions made. In S^] we present the simu- 
lation results, focusing on the general temporal evolu- 
tion from the start until >600ms, a time by which the 
neutrino-driven wind has reached a steady state. We 
then analyse in more detail the various components of 
these simulations. In f^l we discuss the properties of the 
nascent neutron stars, with special attention paid to the 
geometry of the neutrinosphercs. In fj^l we discuss the 
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neutrino signatures, both in terms of luminosity and en- 
ergy distribution. In |3 we focus on the residual material 
lying at low latitudes, forming a quasi-Keplerian disk. In 
S|S| we turn to the energetics of the explosion and describe 
in detail the main component of the simulations at late 
times, i.e., the neutrino-driven wind. In ^ we analyse 
the electron fraction of the ejected material and address 
the relevance of the AIC of white dwarfs for neutron-rich 
element pollution of the interstellar medium. In i jlUI wc 
present the gravitational-wave signal predicted for the 
AIC of our white dwarf models. In Hill we wrap up with 
a discussion of the main results of this investigation and 
present our conclusions. 

2. INITIAL MODELS 

In this section, we present the properties of the AIC 
progenitors selected in our study and summarize the pre- 
sentation in §2 of YL05. 

The general assumption for the construction of 2D pro- 
genitor models for the AIC of white dwarf (YL05) is that 
the resulting structure of the object is essentially inde- 
pendent of its evolutionary history, the only factors that 
matter being the given (final) mass, angular momentum, 
and central density, pc- Additionally, the angular veloc- 
ity distribution ri(r, z), where r is the cylindrical radius 
and z is the distance to the equator, is determined self- 
consistently, given a number of properties identified in 
ID models: 1) the role played by the dynamical shear 
instability, 2) the compression (and spin up) of the sur- 
face layers due to mass accretion which puts its peak 
angular velocity interior to the surface radius, 3) the 
adopted surface rotational velocity value (its fraction of 
the local Keplerian value), and 4) the geometry of the 
angular velocity profile, which we assume to be constant 
on cylinders, i.e., 17 = ^{r). Together with the pres- 
sure/density dependence P — P{p), such rotating stars 
are called barotropic. Note that the criterion used in 
YL05 for the shear rate for the onset of the dynamical 
shear instability is determined using the EOS of Blin- 
nikov et al. (1996). 

The 2D rotating models then correspond to equilib- 
rium configurations iteratively found from trial den- 
sity and angular velocity distributions, using the 
Self-Consistent-Field method (Ostriker & Mark 1968; 
Hachisu 1986), under the constraint that the density p is 
solely a function of the effective potential ^'(r, z), given 
as the sum of the gravitational potential, i.e., 

and the centrifugal potential, i.e., 

e(r) = - / n^{r')r'dr\ 



where R' = \/r'2 + z'2 (Tassoul 2000, YL05). 

In this paper, we select two progenitors with masses 
of 1.46 M© and 1.92 M© and present their global charac- 
teristics in Tabled Both models have an initial central 
density pc equal to SxlO-^^gcm""^. The 1.46-M0 model 
serves as a reference for an object with a moderate initial 
rotational energy T relative to gravitational energy \ W\, 
i.e., r/|iy| — 0.0076, and, indeed, shows only a modest 
initial departure from spherical symmetry, with a polar 



TABLE 1 

Properties of selected AIC progenitors 



M _Rp fleq J T \W\ T/\W\ T/\W 

Mq km km erg-s erg erg initial final 
(1050) (10^0) (lO™) 



1.46 800 1130 0.160 0.7 91.97 0.0076 0.059 
1.92 660 2350 1.092 10.57 126.9 0.0833 0.262 



(i?p) to equatorial (i?oq) radius ratio Rp/Rcq — 0.7. At 
the other end of the white dwarf mass spectrum, the 1.92- 
Mq model has considerable rotational energy, with initial 
T/|W| = 0.0833, and the morphology of the star departs 
strongly from spherical symmetry, with Rp/ R^q = 0.28. 
For later reference, we also provide in Tabled the ratio 
T/|VF| at the end of the simulations. In practice, the in- 
fall of the ambient material causes numerical difficulties 
soon after the start of the simulation. These difficulties 
were resolved by trimming the outer and low-density lay- 
ers of both white dwarf progenitors. While the polar ra- 
dius is hardly affected, the equatorial radius is reduced to 
980 km (down from 1130 km) for the 1.46-M0 model and 
1860 km (down from 2350 km) for the 1.92-Mq model. 
The progenitor mass is, however, reduced by less than 
one part in 10^. Hence, we do not expect any percepti- 
ble effect on the results. 

3. VULCAN/2D SIMULATION CODE 

The simulations discussed in this paper were performed 
with the Newtonian hydrodynamic code VULCAN/2D 
(Livne 1993), supplemented with an algorithm for neu- 
trino transport as described in Livne et al. (2004) and 
Walder et al. (2005). The version of the code used here 
is the same as that discussed in Dessart et al. (2005) 
and Burrows et al. (2005), and uses the 2D Multi- 
Group Flux-Limited Diffusion (MGFLD) method to han- 
dle neutrino transport (see Appendix A of Dessart et 
al. 2005). The MGFLD variant of VULCAN/2D is 
much faster than the more accurate, but considerably 
more costly, multi-angle S'„ variant. Doppler velocity- 
dependent terms are not included in the transport, al- 
though advection terms are. Frequency redistribution 
due to the subdominant process of neutrino-electron 
scattering is neglected. Our calculations include 16 en- 
ergy groups logarithmically distributed in energy from 
2.5 to 220 MeV, take into account the electron and anti- 
electron neutrinos, and bundle the four additional neu- 
trino and anti-neutrino flavors into a "f^" component. 
VULCAN/2D uses a hybrid grid, switching from Carte- 
sian in the inner 20 km to spherical-polar further out. 
In the simulation of the 1.46-M0 model, mapped onto a 
180° wedge, nearly perfect top-bottom symmetry about 
the equatorial plane was maintained during the pre- and 
post-bounce evolution. Thus, for the 1.92-M0 , we lim- 
ited the computational domain to just one hemisphere. 
The thus reduced number of zones was used to increase 
the resolution. To summarize, the 1.46-M0 model uses a 
grid with a maximum resolution in the Cartesian inner 
region of 0.56 km, and a minimum resolution of 150 km at 
a maximum radius of 5000 km, with 121 regularly spaced 
angular zones to cover 180°. The 1.92-M0 model uses a 
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Fig. 1. — Top: Color map of the density, p{r,z) (log-scale), at the start of the simulations for the 1.46-Mq model (left) and the 
1.92-Mq model (right). We overplot the equipotentials corresponding to this equilibrium configuration, whereby the density p{r, z) is 

solely a function of the effective potential ^{r,z) given as the sum of the gravitational potential, i.e., ^{r,z) = —G f ^^^^i^ d?R' , and 
the centrifugal potential, i.e., B(r) = — J Q^{r')r'dr' (see text and YL05 for details). Bottom: Same as above, but this time for the 
angular velocity, fl{r,z). We overplot contours of the temperature T{r,z), chosen to follow the relation T{r,z) = Tc{pc / p{r, z))^'^^ , where 
Tc = 10l° K (Tc = 1.3 X IQl" K) for the 1.46-M0 (1.92-M0 ) model. Both models have an initial central density pc equal to SxlQl" gcm-3. 
Contours shown are linearly distributed between 4x10^ K and lO'^'^ K. 



grid with a maximum resolution in tlic Cartesian inner 
region of 0.48 km, and a minimum resolution of 100 km at 
a maximum radius of 4000 km, with 71 regularly-spaced 
angular zones covering 90°. 

A tricky part of the set up was to choose the properties 
of the "ambient" medium surrounding the AIC model. 
This need arises because our inputs are Lagrangean in 
spirit, while VULCAN/2D employs a Eulerian grid. Ide- 
ally, one would like to have material that merely occu- 
pies the space that will soon, after bounce, be replaced 
by the ejected material following the explosion. A key 
requirement is, thus, that this material have a very low 
pressure, to influence as little as possible the properties 
of the blast, but to allow for a smooth transition from cir- 
cumstellar to ejected material in a given region of the Eu- 
lerian grid. To achieve this, we extended our SHEN EOS 
(Shen et al. 1998) down to very low densities (lOgcm"'^) 



and low temperatures (10^ K), conditions in which the 
medium is actually radiation-dominated and, thus, has 
a pressure that depends mostly on temperature. A suc- 
cessful choice was to adopt a density of 1000 gcm~^ and 
a temperature of 4xl0*K for the ambient medium sur- 
rounding both white dwarf models. 

A major deficiency of the white dwarf progenitor mod- 
els used here is their unknown initial thermal structure, 
which YL05 did not provide. Given the additional dif- 
ficulty in handling low temperatures and high densities, 
we resorted to using a parameterized function of the local 
density, i.e., 

with Tc = 10i°K for the 1.46-Mo model, and Tc = 
1.3 X 10^" K for the 1.92-Mq model. Note that, similarly, 
Woosley & Baron (1992) were forced to set a central tem- 
perature of 1.2x10^"K at the start of their simulation. 
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Figure^recapitulates the basic properties of the white 
dwarf progenitors (left column: 1.46-M0 model; right 
column: 1.92-Mq model) as mapped onto our Eulerian 
grid. In the top panel, we show a color map of the density 
on which we superpose line contours of the effective po- 
tential 5* (defined above), as computed by VULCAN/2D. 
Our computation of the gravitational potential is based 
on a multipole expansion in spherical harmonics up to 
I = 33. We reproduce the fundamental property of these 
fast-rotating white dwarf progenitors, in that the iso- 
pressure surfaces and the equipotentials coincide. In the 
bottom panel, we plot the initial angular velocity field. 
To avoid the distortion of streamlines of the infalling am- 
bient material, the angular velocity is set to zero outside 
the WD progenitor. Note how, from zero, the angular 
velocity rises to a maximum value near the outer equa- 
torial radius in the 1.46-M0 model (left column), while 
it is much higher in the center of the white dwarf in the 
1.92-M0 model (right column). As we will show below, 
rotation has a much bigger impact on the collapse phase 
for the latter configuration. We also overplot line con- 
tours of the temperature for each model, following the 
prescription outlined in the above paragraph. 

Finally, we adopt an initial electron fraction of 0.5; 
the high initial central density permits fast electron cap- 
ture which soon decreases the Yc in the core, leading to 
bounce on a timescale ten times shorter than typically 
experienced in the collapse of the core of massive star 
progenitors (Woosley & Weaver 1995; Heger et al. 2000; 
Woosley et al. 2002). 

4. SIMULATION RESULTS 

In this section, we discuss the general properties of the 
pre- and post-collapse phases for both models at the same 
time. We follow the post-bounce evolution of the 1.46- 
M© (1.92-M0 ) model for 550ms (780ms), for a total of 
520 000 (820 000) timesteps. 

In Figs. 12113 we provide entropy (saturated at 
20kB/baryon, where ke is Boltzmann's constant) color 
maps showing the key events in the evolution of both 
white dwarf models, with starting conditions discussed 
in the previous section and displayed in Fig.^ We com- 
plement these figures with Fig.^for the electron fraction 
evolution at three reference times (left column: 1.46- 
M© model; right column: 1.92-M0 model). Note also 
that we overplot most figures with black arrows repre- 
senting velocity vectors, whose maximum length is set 
to 10% of the width (or the height) of the image. The 
corresponding maximum velocity is then given for each 
image in the figure captions - we also mention if we sat- 
urate the vector lengths. Having a large central density 
of 5xl0^°gcm~^, both models achieve nuclear densities 
after the same time of ~37ms. Differences in the bounce 
properties are attributable to the initial inner angular 
velocity distributions (Fig. Compared to the 1.46- 
Mq model, the faster rotator has a lower maximum den- 
sity at bounce (slightly shifted from the grid center to 
r ~lkm), i.e., 2.2 instead of S.lxlO^^'gcm"'^. It mani- 
fests an oblate, rather than a spherical, inner region (the 
inner tens of kilometers) of low entropy (^Ifce/baryon; 
visible in red). The deleptonized material in this inner re- 
gion is, however, aspherically distributed in both models, 
although more so in the 1.92-M0 model, with the lower 
Yc material lying in a disk structure along the equatorial 



direction. The flatter density gradient in this direction 
and the longer dwell time near the neutrinosphere con- 
spire to produce this enhanced deleptonization. At core 
bounce, the low-density outer parts of the white dwarf 
progenitor have not yet started to infall and the progen- 
itor still retains its original shape. The collapse of the 
inner regions, however, creates a rarefaction wave that 
triggers the infall of the outer material, with a magni- 
tude that is more pronounced along the poles due to the 
compact structure of the white dwarf in these directions. 

Although we do not observe a prompt explosion (i.e. 
occuring on a dynamical timescale of just a few mil- 
liseconds) , the shock progresses slowly outwards without 
conspicuously stalling. This is slightly different from the 
core-collapse simulations recently published, whereby the 
shock systematically stalls (Burrows et al. 2005; Buras 
et al. 2005ab). In the equatorial direction, facilitated 
by the centrifugal support of the infalling material, the 
shock progresses outwards steadily, and is faster in the 
faster rotating model, reaching a few hundred kilometers 
after ^lOOms. Given the initial constant rotation rate 
on cylinders (see |5Jj ^^e radial inflow of mass brings an- 
gular momentum to the equatorial (z 0) region, which 
becomes more strongly supported, facilitating further the 
progress of the shock at low latitudes. This exacerbates 
the non-spherical development of the shock structure in 
both models. In its wake, we see a few large-scale whirls, 
resulting from the generation, at mid-latitudes and due 
to shock passage, of large vortical motions. 

In the polar direction, centrifugal support is absent, 
but the stalling of the shock is prevented by the quickly 
decreasing accretion rate established by the steeper den- 
sity gradient, reduced polar radius, and smaller mass 
budget. Hence, although the shock slowly migrates along 
the equatorial regions, it soon traverses the surface layer 
of the white dwarf along the poles and escapes outwards 
into the ambient medium. Due to the strong aspheric- 
ity of the progenitor, this occurs only 70 ms after core 
bounce in the 1.92-M0 model, ~30ms earlier than for 
the 1.46-M0 model. The ambient medium is then swept 
up by this blast, whose opening angle is constrained by 
that of the uncoUapsed disk of the progenitor. The out- 
flow expansion rate is larger 30-40 degrees away from 
the poles than right along the pole, coinciding in Fig. 0] 
with the lower Y^ material. The effect is large for the 
1.92-M0 model, giving a butterfly shape in cross section 
to the faster expanding portions of the outflow. Off-axis 
material has more rotational kinetic energy available to 
convert to 2D planar (r, z) kinetic energy as it streams 
outward, reducing its rotational velocity while preserv- 
ing angular momentum, thereby resulting in an enhanced 
acceleration compared to that of the material situated 
along the poles and lacking rotational energy. 

As the shock expands, it wraps around the disk, typ- 
ically with a speed near that of the local sound speed. 
Nearer the pole, the outflow sweeps along the pole-facing 
side of the pole-excavated white dwarf, entraining surface 
material and effectively loading the outflow with more 
mass, causing unsteady fallback onto the neutron star. 
For the 1.46-M0 progenitor, the shock completely wraps 
around the low-latitude regions of the white dwarf, and 
finally emerges from the outer equatorial regions, as wit- 
nessed by an outward- moving entropy jump. By 250 ms 
after bounce, the shock reaches a few thousand kilome- 
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Fig. 2. — Slide sequence of color maps of the entropy covering the post-bounce evolution of the 1.46-Mq progenitor, with velocity vectors 
overplotted as black or white arrows. The entropy, reaching up to SOfca/baryon at later times, has been saturated at 20fcB/baryon in all 
panels to enhance the contrast. From top to bottom, left to right, the vector of maximum length corresponds to a velocity magnitude of 
52500 km s-l(infall), 28500 km s"! (infall), 24000 km s^^ (infall), 20500 km s'^ (outflow), 16500 km s~l (outflow), and 34000 km s"! (outflow). 
We indicate the time after bounce in the bottom-left corner of each panel. Note also that the spatial scale varies between panels, with 
widths in the range 300 km to 10000 km. See text for details. 
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Fig. 3. — Slide sequence of color maps of the entropy covering the post-bounce evolution of the 1.92-Mq progenitor, with velocity vectors 
overplotted as black or white arrows. The entropy, reaching up to 50 kg /baryon at late times, has been saturated at 20 kg /baryon to enhance 
the constrast. From top to bottom, left to right, the vector of maximum length corresponds to a velocity mag nitude of 53000 km s"! (infaU), 
36500 kms~l(infall), 29000 km s'^ (infall), 32000 km s"^ (outflow), 32500 km s~^ (outflow), and 34500 km s"^ (outflow). We indicate the 
time after bounce in the bottom-left corner of each panel. Note also that the spatial scale varies between panels, with widths in the range 
300 km to 8000 km. See text for details. 
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Fig. 4. — Slide sequence of color maps of the electron fraction, 1^, covering the post-bounce evolution of the 1.46-Mq (left column) 
and 1.92-Mq (right column) progenitors, overplotted with velocity vectors, as black arrows. We also superpose isodensity contours as 
white lines, spaced every decade, starting at lO^^gcm""^ (top row) or lO^^gcm"^ (bottom two rows). From top to bottom, left to right, 
the vector of maximum length corresponds to a velocity magnitude of 43500 km s~^ (infall), 36500 km s~^ (infall), 25500 km s~^ (outflow), 
19500 km s~^ (outflow), 34000 km s~^ (outflow), and 34500 km s~^ (outflow). We indicate the time after bounce in the bottom-left corner 
of each panel. The spatial scale varies between panels, with widths in the range 500 km to 10000 km (8000 km) for the 1.46-M0 (1.92-M0 ) 
model. See text for details. 
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Fig. 5. — Top: Time evolution of the polar density profile for the 1.46-M0 (left column) and 1.92-Mq (right column) progenitor models. 
Bottom: Same as above, but for the equatorial direction. Note the steeper density profiles at all times along the pole compared to the 
equatorial direction, highlighting the contrast between the wind (polar direction) and disk (equatorial direction; within a 2000 km radius 
at 770 ms for the 1.92-Mq model). See text for discussion. 
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Fig. 6. — Same as for Fig. |21 but this time for the radial velocity Vr. Note the large contrast between the polar and equatorial 
pre-bounce (infall) velocities for the 1.92-M0 model, the fast neutrino-driven wind at late times along the poles, overtaking the slower 
expansion of the shock generated at core bounce. Along the equator, material shows slow (1.46-M0 model) or absent (within a 2000 km 
radius, 1.92-M0 model) expansion, belonging to a quasi-Keplerian disk (JJJ. See text for discussion. 



ters, assumes a near spherical shape, and the entire white 
dwarf material outside the newly-formed neutron star 
flows nearly radially outward. In the high- rotation (1.92- 
Mq ) progenitor, the white dwarf possesses a lot more 
mass at near zero-latitude and this confines more drasti- 
cally the emerging shock along the poles. As the shock 
migrates outwards, it opens up; it does wrap around the 
progenitor, but much later than when the shock escapes 
in the polar directions. Despite the reduced ram pres- 
sure associated with centrifugal support, the shock stalls 
a few hundred milliseconds after bounce along all near- 
equatorial directions. 

Within 100-200 ms, the newly-formed neutron star has 
a mass of ~1.4Mq , similar in both models despite the 
0.5 Mq difference in progenitor mass. Note that the 
large neutron star asphericity and the sizably lower den- 
sity at the neutrinosphere for near-zero latitudes makes 
this mass definition ambiguous at such early times, es- 
pecially in the 1.92-M0 model. Indeed, the neutron 
star is not clearly distinguishable from the surround- 
ing equatorial material, so imposing either a density 
cut of 10^°~^^gcm~^ or a radius cut in defining the 
newly-born neutron star appears arbitrary when deter- 
mining the residual mass. In the 1.46-M0 model, about 
0.06 Mq remains outside of the neutron star, mostly 
in the equatorial disk region; in the 1.92-M0 model, 
0.6 Mq is now lying in this disk-like structure. The rest 
of the initial mass is outfiowing material, which, if se- 
lected according to an outward radial velocity discrimi- 



nant of 10000 kms ^(comparable to the escape velocity 
at 3000 km), reaches 4xlO-3MQfor the 1.46-Mq model 
and 3x 10-3 Mq for the 1.92-Mq model. These various 
components are documented in more detail in §|3-|51 

The late-time evolution of both models is character- 
ized by a strong neutrino-driven wind that sets in about 
300 ms after bounce, replenishing the grid with denser 
material (on average lO^gcm^^) and large velocities 
(with a maximum of 30000 km s^^ along the poles). The 
properties of the neutrino-driven wind are very angle- 
dependent, the density changing by 30% between the 
pole and the 40° latitude, while the radial outflow ve- 
locity varies by a factor of 3 in the 1.92-Mq model. The 
latitudinal dependence of the mass flux per unit solid an- 
gle is therefore dominated by a variation in asymptotic 
velocity. We will discuss this result in more detail in 
SjSl By the time we stop the simulations, at 550 ms and 
780 ms for the 1.46-Mq and 1.92-Mq models, all the am- 
bient medium originally placed around the white dwarf 
progenitor has been swept away by the neutrino-driven 
wind, which occupies all the space outside the neutron 
star and the disk. The electron fraction of the material 
ejected in the original blast is close to 0.45-0.5, while 
subsequently, in the neutrino-driven wind, the values are 
lower, with a pronounced decrease towards lower lati- 
tudes (note, however, the high Ye right along the pole; 
see fQ- 

In Fig. [SI we recapitulate for both models the evolu- 
tion described above by showing equatorial and polar 
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slices of the density as a function of time. Striking fea- 
tures are the distinct polar and equatorial surface radii, 
the fast infall of the inner regions to nuclear densities, 
the slow plowing of the shock along the equatorial di- 
rections, superseded in radial extent and velocity by the 
shock in the polar direction as the surface mass shells col- 
lapse in, and finally the emergence of a sequence of den- 
sity kinks associated with the birth of the fast neutrino- 
driven wind that sweeps away the previously shocked ma- 
terial that did not leave the grid. Similarly, in Fig. 1^1 
we show slices of the radial velocity, Vr, along the po- 
lar (top row) and equatorial (bottom row) directions for 
the 1.46-Mq (left column) and 1.92-Mq (right column) 
models. Notice the much larger infall velocities, similar 
along the poles and the equator for the 1.46-M0 model, 
but with a strong latitudinal dependence in the 1.92- 
M0 model. In that model, the speed contrast between 
the polar and equatorial directions is ^30000 kms~^. 
Overall, the evolution is more rapid along the poles 
than on the equator, with larger asymptotic velocities 
(30000 kms~^compared with 10000 kms~^), and with 
the establishment of a quasi-stationary outflow at late 
times along the pole. These radial slices offer a means 
of better interpreting the fluid velocities, depicted with 
vectors, in most color maps shown in this paper. We also 
show isodensity contours in most color maps to provide 
some feeling for the density distribution. 

Having described the general properties of the two 
simulations of the AIC of a 1.46-M0and 1.92-M0 white 
dwarf, we now address more speciflc issues, covering the 
properties of the nascent neutron star {^^, the neutrino 
signatures (SQ, the properties of the residual disk and 
the angular momentum history (13, the neutrino-driven 
wind and the global energetics {^^, the electron fraction 
of the ejected material {^^, and, finally, the gravitational 
wave signatures fiil()|l. 

5. NEUTRON STAR PROPERTIES 

The white dwarf progenitors discussed in this paper, 
due to their evolution to high central densities, high ro- 
tational kinetic energies, and high mass, are distinctive 
in that their cores always collapse to form neutron stars 
rather than being disrupted by the explosive burning of 
carbon and oxygen. 

We find that neutron stars formed from the AIC of 
the progenitor white dwarfs used in this work are very 
aspherical (see, for example, Walder et al. 2005; Liu & 
Lindblom 2001; Janka & Monchmeyer 1989ab), although 
there are significant differences in evolution after bounce 
between the two models. In Fig.[7| we show color maps of 
the density field 59 ms after bounce (top row) , as well as 
for the last time computed {t — 570 ms after bounce; bot- 
tom row) in the 1.46-M0 model. To render more striking 
the level of asphericity of the neutron star "surface," we 
overplot the neutrinospheres R^{e^,r, z), adopting the 
definition 

/•OO 

r(i?,(e„r,z)) = / K,{p,T,Y,)p{r' , z')dR' = 2/3 , 

where K^{p,T,Ye) is the combined material absorption 
and scattering opacity to neutrinos, and the integration 
is carried out along radial rays, with R' = \/r'^ + z'^, 
using 30 equally spaced latitudinal directions per quad- 
rant. Strictly, this definition is most appropriate for 



the photosphere/neutrinosphere in a plane-parallel at- 
mosphere, but it gives a sense of the asphericity of the 
collapsed core. Note that the above expression contains 
a dependence both on the neutrino flavor and the neu- 
trino energy e^. In the left column of Fig. 13 line con- 
tours correspond to such neutrinosphere radii as a func- 
tion of energy group, bracketing the peak of the neutrino 
energy distribution at the neutrinosphere, i.e., between 
2.5 and 46MeV. Material opacity to neutrinos increases 
with the square of the energy, so that higher-energy 
neutrinos have larger neutrinospheres. Here, for the 
1.46-M0 model, these radii vary from ~30 to ^120 km 
along the equatorial direction, with little departure from 
sphericity (30% lower values are obtained in the polar 
direction). In the right panel, line contours correspond 
to the neutrinosphere for the three different flavors at 
a neutrino energy (e^) of 12.5 MeV, associated with the 
peak of the energy distribution at inflnity. Note that 
this is approximate, since the neutrino energy distribu- 
tion hardens with time and manifests a latitudinal depen- 
dence (see As in standard ID and 2D core-collapse 
computations, we find that the electron neutrinos de- 
couple from matter at larger radii than the Ve and 'V;^" 
neutrinos. Here, the former decouple at 90 km (60 km) 
along the equator (pole) , the latter two further in but at 
a similar radius of 70 km (50 km). The neutrinospheres 
show a similar shape for all three flavors (and all energy 
groups), reflecting the corresponding asphericity in the 
density fleld. 

In the bottom-row panels, we reproduce the above for 
the last time in the 1.46-M0 simulation. The departure 
from sphericity is now considerable, with both an oblate- 
ness and a strong pinching of the neutrinospheres along 
the polar directions. Along the equatorial direction, the 
radial spacing between neutrinospheres of consecutive 
and higher energy groups has increased, and the lower 
(higher) energy groups decouple further in (out) than in 
the previous snapshot, with neutrinospheres between 20 
and 150km (from 2.5 to 50 MeV), 80 and 50km (for the 
Ve neutrino, and v^j^^v^^^ neutrinos, respectively). Along 
the polar direction, neutrinospheres of all neutrino en- 
ergy groups (and all neutrino flavors at 12.5 MeV) shown 
reside in a narrow range of radii between 20 and 30 km 
(22 to 25 km). Here again, the neutrinospheres depicted 
follow very closely the contours of density, which is the 
primary factor controlling the neutrino optical depth. 

In Fig.|Sl we duplicate Fig. [T] (note the different spatial 
scale) for the 1.92-M0 model, showing the same quan- 
tities both early after bounce (59 ms) and significantly 
later at 775.5 ms after bounce. There are numerous 
differences with the 1.46-M0 model. First, the neutri- 
nospheres are aspherical even right after bounce (top 
row), with equatorial (polar) radii larger (smaller) by a 
factor of 2-3 compared with those in the 1.46-M0 model. 
All flavors reveal similar neutrinosphere locations. 

In the 1.46-M0 model, the later ratio of the equatorial 
and polar radii is 2.5:1, irrespective of energy group and 
flavor. This becomes 15:1 in the 1.92-M0 model, and is 
thus a considerable departure from sphericity; the faster 
rotating model has a neutrinosphere radius of just 14 km 
along the pole, but 215 km along the equator. This is the 
most conspicuous difference with the essentially spheri- 
cal neutron stars seen in non-rotating simulations of the 
more standard core collapse of massive stars (Keil et al. 
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Fig. 7. — Top: Color map of the density at 59 ms after bounce for the 1.46-Mq model. Left: Neutrinosphere radii for the neutrino for 
ten energy groups around the peak of the emergent energy distribution, i.e., at 2.5, 5.5, 7.5, 9.5, 12.5, 16, 21, 27, 36, and 46MeV. Right: 
Neutrinosphere radii at 12.5 MeV for the three neutrino flavors: Uc (solid lines), Uc (dotted line), and 'V^" (dash-dotted line, which overlaps 
with the dotted line). Bottom: Same as above, but for the last computed time at 570 ms after bounce. The vector with the maximum 
length corresponds to a velocity magnitude of 46300 km s^^ (infall; top panels) and 9440 km s~^ (outflow; bottom panels). The fixed spatial 
coverage of the display is 200 x 200 km^ . 



1996; Swesty & Myra 2005; Buras et al. 2005a; Dessart 
et al. 2005), with neutrinosphere radii of the order of 
20-30 km at comparable times after core bounce. 

The very aspherical neutron stars formed through the 
AIC of a white dwarf make the determination of the 
neutron star mass somewhat ambiguous. Rather than 
taking the enclosed mass within a given spherical ra- 
dius, we compute the total mass from all regions above 
a given mass density. For a density cut of lO^^gcm"^, 
we obtain neutron star masses of 1.42 Motor the 1.46- 
M© model, and 1.5 Motor the 1.92-M0 model. However, 



if we adopt a density cut of lO'^^gcm"'^, the neutron 
star masses are, respectively, 1.39 M© and 1.30 M© . The 
higher-mass progenitor model has now a smaller neutron 
star mass, reflecting the strong asphericity of the den- 
sity field. While we might associate the neutrinosphere 
with the neutron star surface and with a standard mass 
density of lO'^^ gcm~'^, such an association in the present 
fast rotating neutron star is inappropriate, since the neu- 
trinospheres extend well into regions where the density 
is <10^"gcm~^. These neutron star mass values are 
reached at 100 ms after core bounce and remain essen- 
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Fig. 8. — Top: Color map of the density at 59 ms after bounce for the 1.92-Mq model. Left: Neutrinosphere radii for the Uc neutrino for 
ten energy groups around the peak of the emergent energy distribution, i.e., at 2.5, 5.5, 7.5, 9.5, 12.5, 16, 21, 27, 36, and 46MeV. Right: 
Neutrinosphere radii at 12.5 MeV for the three neutrino flavors: Uc (solid lines), (dotted line), and 'V^" (dash-dotted line, which overlaps 
with the dotted line). Bottom: Same as above, but for the last computed time at 775.5 ms after bounce. The vector of maximum length 
corresponds to a velocity magnitude of 36500 kms~^(infall; top panels) and 15500 km s^^ (outflow; bottom panels). The fixed spatial 
coverage of the display is 500 X 500 km^ . 



tially constant. The neutrino-driven wind that appears 
after a few 100 ms decreases the neutron star mass at a 
rate of just a few 10"'^ M0 s~^ (see ISJ. Note also that 
the enhanced centrifugal support in the faster-rotating, 
higher-mass model leads to bounce at a 30% lower max- 
imum density compared with the 1.46-Mq model, and 
both a reduced and a delayed mass accretion rate along 
the equatorial direction compared with what would pre- 



vail in the absence of rotation. At the end of each sim- 
ulation, the total neutron star angular momentum is 
1.35X lO''^ erg-s (1.13x10''^ erg-s) for the 1.46-M0 model, 
using the density cut at 10^°gcm~^ (lO^^gcm"'^) , and 
4.57x10*9 erg-s (2.79xl0''9 erg-s) for the 1.92-M0modeL 
However, the accretion rate at a given Eulerian radius 
is higher and longer-lived at smaller latitudes, because 
of the larger amount of mass available and the flatter 
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density profile in those regions. Overall, the presence 
of a massive accretion disk in the fast rotating model 
complicates the definition of the neutron star mass at 
such early times. Evolution over minutes/hours/days 
will likely lead to significant accretion onto the neutron 
star, resulting in a much higher final mass. 

The final rotational to gravitational energy ratio 
T/\W\ is 0.059 for the 1.46-Mo model and 0.262 for the 
1.92-Mq model. These values are large and for the lat- 
ter model, large enough to cause the growth of secular 
and perhaps even dynamical instabilities (Tassoul 2000). 
Using realistic post-bounce configurations for a rotating 
massive star progenitor, Ott et al. (2005a) find a dynam- 
ically unstable spiral mode for T/|T4^| as low as ^0.08. 
Thus, it is likely that the PNS structures found here, es- 
pecially for the 1.92-Mq model, would develop some non- 
axisymmetric instability that would cause, among other 
things, outward angular momentum transport. 

These results are significantly different from those of 
Fryer et al. (1999), who obtained a successful explosion 
<100ms after bounce, a PNS mass of ~1.2M0 , and a 
~1 s period at ~200ms after core bounce. Such different 
conclusions stem from their adoption of slow, solid-body 
progenitor rotation, with most of the angular momen- 
tum stored in the outer mantle and blown away by the 
explosion, rather than being accreted by the PNS. 

6. NEUTRINO SIGNATURES 

The first observational signature of an AIC explosion 
would be the copious emission of neutrinos immediately 
after core bounce. As discussed above, the properties 
of the bounce of the core and the nascent neutron star 
are close enough to those obtained in simulations of the 
core collapse of massive progenitors that one expects a 
neutrino signal with a somewhat similar evolution and 
character (see, for example, the predictions for the 11- 
Mq model of Woosley & Weaver 1995 in Dessart et al. 
2005). 

In Fig. 13 we show the neutrino luminosities (on a log 
scale) for the 1.46-M0 (left panel) and 1.92-Mo (right 
panel) models, with distinct curves for the different neu- 
trino flavors (solid line: Ve] dashed line: Ve] dash-dotted 
line: "f/j")) as well as different colors for the equatorial 
(black) and polar (red) directions. These are luminosi- 
ties in the sense that the flux in each direction is scaled 
by 47ri?^ where i? is a spherical radius, of 250 km for 
the 1.46-M0 model and 400 km for the 1.92-Mq model 
(chosen to be well above the neutrinosphere for the en- 
ergy at the peak of the neutrino distribution). They 
correspond to the total luminosity that would have ob- 
tained had the selected directional flux been the same 
in all directions. The temporal evolution of the various 
fluxes for both models is comparable. The total neu- 
trino luminosity reaches a maximum of 5.2x10^^^ ergs"^, 
mostly due to the Ve neutrino contribution, and decreases 
to ~4x 10^^ ergs~^ at 500ms after bounce in the 1.46- 
Mq model, with a further 30% decrease for the 1.92- 
M© model. At later times, the main reason for this dif- 
ference is the much lower 'V^" neutrino luminosity in 
the 1.92-Mo model. This reduction has been seen and 
discussed by Fryer & Heger (2000) in the context of the 
collapse of rotating cores of massive progenitors. The 
smaller core densities (weaker bounce) achieved in mod- 
els with fast rotation lead to smaller temperatures and, 



consequently, smaller neutrino emission, with a larger ef- 
fect for the and Vr neutrinos (grouped under the name 
"f^" here). So, while the globally lower neutron star 
densities in the fast-rotating model induce a reduction 
in neutrino luminosity compared to the 1.46-M0 model, 
the same effect introduces a latitudinal variation of neu- 
trino fluxes in the faster rotating model, with fluxes, ir- 
respective of neutrino flavor, larger by a factor of about 
two along the pole than along the equator (note that 
at a radius of 250 km, the difference is higher and on 
the order of three). This variation, not discussed by 
Fryer & Heger (2000), results from the further varia- 
tion of the neutrinosphere temperatures with latitude 
within a given rotating model. In both models, but more 
so in the 1.92-Mq model, the temperature gradient and 
the temperatures are reduced at the neutrinosphere (for 
a given energy group) along the equator compared to 
the poles, irrespective of energy group and flavor, as is 
clearly visible in Fig. E| For example, we see that the 
temperature on the 10^°gcm~'^ contour is 4MeV in the 
polar direction and 0.5 MeV in the equatorial direction. 
Within the neutrinosphere region, we find fluid velocities 
that are oriented preferentially in the ^-direction, along 
cylinders, illustrating the weak or absent convection that 
results from the stabilizing specific angular momentum 
profile (Fryer & Heger 2000; Heger et al. 2000; Ott et 
al. 2005b; see also SEj). Along a given angular slice, the 
temperature has a maximum at mid-latitudes, caused by 
enhanced [ve) neutrino energy deposition in this direc- 
tion. These regions offer a tradeoff, since the fiux is still 
higher at such latitudes than along the equator, while 
the density is relatively higher than along the poles (see 
The latitudinal variations seen in the collapsed mod- 
els of AIC progenitors are extreme, and, indeed, for the 
slower rotation rates typically obtained for massive-star 
core-collapse progenitors (Heger et al. 2000), a modest 
anisotropy is found instead (Walder et al. 2005). 

We document further the neutrino signatures of the 
AIC of white dwarfs by showing in the top row of Fig. lllL 
for the 1.46-Mq (left) and 1.92-M0 (right) models, the 
time- integrated neutrino emission at infinity for each fla- 
vor, as a function of neutrino energy. In the bottom three 
panels of each column, we show the individual neutrino 
distributions at three representative times of the simu- 
lations (at bounce, halfway through the simulation, and 
at the last simulated time). The overall flux level and 
hardness of the energy distribution arc higher along the 
polar direction, the variation at a given time towards 
higher latitude mimicking the time evolution seen for 
non-rotating core collapse simulations of massive star 
progenitors. We compute the average neutrino energies, 
here defined as 



J d£^F^{£^, R) 

For the 1.46-M0 model [R — 250km), we obtain similar 
values to within 2-3% along the pole and along the equa- 
tor with {e^^) = 10 MeV, (e^J = 15 MeV, and (e^^J = 
24 MeV. For the 1.92-Mq model (i? = 600 km), we obtain 
systematically lower values than for the 1.46-Mo model, 
and for lower latitudes. Along the equator (pole), we ob- 
tain (e^J = 9MeV (lOMeV), (e^J = 14MeV (16MeV), 
and {eu) = 16 MeV (21 MeV). Note that this definition 
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Fig. 9. — Flavor- and angle-dependent neutrino luminosities for the 1.46-M0 (left) and 1.92-Mq (right) progenitors, for an adopted 
radius of ij = 250 km for the former and R = 400 km for the latter. Luminosities shown as solid, dash-dotted, and dashed lines correspond 
to J^^ dQR^ Fv{R, 9, £,j) for the Ue, Pe, and 'V^" neutrinos, repectively, and in black for the equatorial direction {9 = 7r/2) and in red 
for the polar direction {6 = 0). Of particular relevance to the asymptotic electron fraction of the ejecta is the higher electron-neutrino 
luminosity compared to that of the anti-electron neutrinos, by a factor of at least ten in the 1.92-M0 model (right). Note that the vertical 
scale extends down to 10^^ ergs~^ for the 1.46-Mq model and to 10^" ergs~^ for the 1.92-M0 model. 




Fig. 10. — Color map of the matter temperature for the 1.92-M0 model at the last snapshot. We overplot the contours of the density 
every decade, starting with 10^*gcm~^ for the innermost line. Left: View of the innermost region (~20km), showing the latitudinal 
variation of the temperature even on the innermost density contour. The velocity field reveals the presence of fluid motions along cylinders 
due to rotation (see J7J, thereby preventing convection. A ten-fold magnification of the velocity- vector length, saturated at a value of 
9130 kms~^, has been applied to better render the fluid motions. Right: Zoom out of the left slide, showing the large temperature 
variation along the 10^""^^ gcm""^ density contour, which corresponds reasonably well to the neutrinosphere location at the peak of the 
neutrino energy distribution (see Fig.|SJ. The vector of maximum length corresponds to a velocity magnitude of 11300 kms~^. 



of the neutrino "average" energy is similar to that found 
in Thompson et al. (2003), who used the mean intensity 
in place of the flux F,^. Since we are close to free- 
streaming regimes at the adopted radii for the peak of 
the neutrino energy distribution, the two are equivalent. 

In the 1.46-M0 model, and at late times, we see that 
the neutrino display is more dramatic (total flux is twice 
as high) and is characterised by a harder spectrum than 
in the 1.92-M0 model. This is due to the more com- 
pact and, thus, hotter neutrinospheres of the neutron 
star formed by the lower-mass white dwarf progenitor. 

7. ROTATION AND THE REMNANT DISK 



Let us now turn our discussion to the angular momen- 
tum and angular velocity budget and profiles in our sim- 
ulations. As discussed in §2]|31 we find that the early 
neutron stars have comparable masses in the two simula- 
tions, the rest residing not so much in the outflow than in 
a substantial amount of "circum-neutron star" disk ma- 
terial, rotating fast, but having little outflow or inflow 
velocity (see Fig. El. 

In Fig. El we plot a temporal sequence of the equato- 
rial radial profile of the angular velocity (top row) and 
specific angular momentum (bottom row), for the 1.46- 
Mq model (left column) and 1.92-Mq model (right col- 
umn). In the 1.46-Mq model, the central angular velocity 
is ~0.1rads~^ (or a period P = 63 s) at the start of the 
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Fig. 11. — Top Left: Time integral, until 570 ms after bounce, of the neutrino luminosity, set equal to ^^B? Fv{R,9)dQ, {9 is the polar 
angle, R is spherical radius, and f2 is the solid angle; see JSJ, at H = 250km, for each neutrino flavor (black: Ve\ blue: Ue\ red: 'V^"), 
and along the polar (solid line) and equatorial (dotted line) directions. Bottom left panels: Neutrino energy distribution, for each neutrino 
flavor, at 10ms (black), 280ms (blue), and 560ms (red) post-bounce times, for both equatorial (dotted) and polar (solid) directions. Right 
column: Same as left column, but this time for the 1.92-M0 model, using a radius R = 600 km. In the bottom three panels, the times 
selected are now 10ms (black), 380ms (blue), and 760ms (red) after bounce. Note the much larger contrast between the equatorial and 
polar fluxes and energy distributions. 
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Fig. 12. — Top: Time evolution of the equatorial angular velocity profile for the 1.46-M0 (left) and 1.92-Mq (right) progenitor models. 
Bottom: Same as above, but this time for the specific angular momentum j. Note its outward increase, which acts as a stabilizing infiuence 
on convective instabilities (Heger et al. 2000), as well as the rise, from zero, to larger values after ~200ms as the neutrino-driven wind wraps 
around the white dwarf and brings rotating material to low latitudes. For all panels, we overplot as a dotted black line the corresponding 
Keplerian value, i.e., f2KGplorian(^) = \J GM j r'^ . See text for details. 
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Fig. 13. — Color map of the logarithm of the angular velocity, at the last simulated time (570 ms for the 1.46-M0 model, and 775 ms 
for the 1.92-Mq model), out to a maximum radius of 2000 km, for the 1.46-Mq (left) and 1.92-Mq (right) progenitors. The initial angular 
velocity profiles f2(r, z), functions of the cylindrical radial coordinate only (YL05 ), a re given for the equatorial direction in Fig. 1121 (bottom 
curve). To enhance the visibility, we have reduced the colorbar range (see Fig. 1121 for the unsaturated values). We also overplot in white 
the corresponding contours of the density for every decade between 10** and 10^" gcm~'^. The vector of maximum length corresponds to a 
velocity magnitude of 29000 km s~^ (outflow) in both panels. The spatial scale is 4000 km in width. 



simulation (initial conditions in the YL05 progenitor), 
and ^lOOOrads"^ (P = 6.3 ms) at the end the simula- 
tion, a spin-up factor of 10000. In the 1.92-M0 model, 
we start with a much higher angular rotation rate of 
'--^20rads~^ (P = 0.3 s), but the final values are compara- 
ble with those of the 1.46-M0 model, being ~2800rads~^ 
(P = 2.2 ms). Thus, both simulations lead to the forma- 
tion of a neutron star with a period of a few millisec- 
onds, although we expect the neutron star formed in the 
1.92-M0 model to further accrete mass and angular mo- 
mentum, which may spin-up the residue to even shorter 
periods. The general angular velocity and specific angu- 
lar momentum profiles for both models are quite similar. 
Despite wiggles observed in the 1.46-M0 model in the in- 
ner 10 kilometers (which we associate with slight numer- 
ical artifacts along the axis - this problem is not present 
in the 1.92-M0 model, whose grid covers only 90°), the 
neutron star is close to solid-body rotation out to 30 km, 
showing a steady and smooth decline with radius beyond. 
In all four panels, we overplot as a broken and black 
line the corresponding local Keplerian angular velocity, 
f^Kcpicrian(?') — \/GM/t^ , whcrc G is the gravitational 
constant and M is the mass interior to the radius (cylin- 
drical, or spherical). (For these plots, we employ the 
corresponding neutron star mass.) The angular veloc- 
ity or specific angular momentum profiles beyond 30 km 
graze the corresponding line for the Keplerian value, be- 



ing always lower by a few tens of percent. In the 1.46- 
M0 model, the profiles evolve significantly toward this 
Keplerian limit, angular momentum being gained along 
the equatorial direction through the radial infall of non- 
zero latitude material. The constancy of the angular ve- 
locity with z allows a significant gain from such infall. 
In the 1.92-M0 model, the rotational properties along the 
equator are originally closer to the Keplerian values, but, 
accordingly, evolve little. In both models, angular mo- 
mentum is transported outwards, first in material blown 
away by the shock wave initiated at core bounce, and 
then in the neutrino-driven wind. This occurs as the 
outflowing material wraps around the progenitor white 
dwarf and eventually meets along the equator. There is 
no "physical" viscosity in the code that would permit a 
proper modeling of the accretion disk. Mass accretion 
should occur in partnership with outward transport of 
angular momentum over a longer timescale, yet to be 
determined. 

In the 1.92-M0 model, at the end of the simula- 
tion, the near-Keplerian disk extends from 30 km out to 
^1800 km, covering a range of densities (temperatures) 



from 10^^ g cm 
0.1 MeV; Vig.^ 



down to 10 gem ^ (3 MeV down to 



8. ENERGETICS AND THE NEUTRINO-DRIVEN WIND 

Given that it leads to the formation of a ^1.4- 
M0 neutron star, the AIC of a 1.4-2.0 M0 white dwarf 
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Fig. 14. — Time evolution of outflow gravitational (blue), thermal (cyan), (r, x)— kinetic energy (red), and rotational (green) energy (on 
a log scale) for the 1.46-Mq model (left) and 1.92-Mq model (right), accounting for all the material with radial outflow velocity higher 
than 10000 kms~^. The explosion energy at the end of the simulation is 2.7xl0*^erg for the 1.46-Mq model and 2xl0'*^erg for the 
1.92-Mq model, with final expected energies of ~10^''erg and ^SxlO^^erg, in the same order. See text for details. 
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Fig. 15. — Top: Time evolution of the total mass loss rate 
through the outer grid radius, dMf^^^^^/dt, in the 1.92-M0 model. 
Bottom: Latitudinal dependence of the radial velocity (solid line) 
and density (broken line) in the 1.92-M0 model at the latest time 
computed, when the neutrino-driven wind has reached a steady- 
state. 



is expected to result, in the case of a successful explo- 
sion, to an outflow of modest mass. Furthermore, due to 
the similarity with the core collapse of massive star pro- 
genitors, their explosion kinetic energy should be lower 
than the --10^^ erg inferred, e.g., for SN1987A (Arnett 
1987). 

As discussed in 21 ^-nd the total sum of the neu- 
tron star and disk masses is very close to the original 
progenitor mass, leaving typically a few 0.001 Mq for the 
ejected material. Integrating all the mass that has left 
the grid over the course of the simulation, as well as all 
the material outflowing with a positive radial velocity 
greater than 10000 km s~^ (which is of the order of the 
escape velocity at a radius of 3000 km) we find a value 
of 4 X 10-3 Mq for the 1.46-M0 model and 3x 10"^ y^^ foj. 
the 1.92-M0 model. In Fig.^1 we show the evolution of 
the corresponding gravitational (blue), thermal (cyan), 
and kinetic (2D planar: red; rotational: green) energies 
for this outflowing material as solid lines, including the 
total energy as a black dotted line. At the last computed 
time, the total energy is indeed lower than that inferred 
for standard core collapse. Adopting a radial velocity 



cut of 10000 kms"^, we find an energy at the last sim- 
ulated time of 2.7xl0''^erg for the 1.46-Mo model and 
2xl0"'^erg for the 1.92-M0 model. Note that energy is 
still being pumped into the wind by the slowly decaying 
neutrino luminosity emanating from the neutron star; 
the trend of the total energy curve suggests that the to- 
tal energy of the explosion will be 2-3 times higher, thus 
^lO^^erg for the 1.46-Mo model and ~5xl0'*^erg for the 
1.92-Mq model. 

This is over one order of magnitude smaller than the 
explosion energy inferred for normal core-collapse su- 
pernovae. The AIC of white dwarfs is likely to lead 
generically to underenergetic explosions because there 
is too little mass to absorb neutrinos, most of it being 
quickly accreted while the rest is centrifugally-supported 
at large radii, far beyond the region where there is a 
positive net gain of electron-neutrino energy. Interest- 
ingly, similar underenergetic explosions are obtained by 
Kitaura et al. (2005) and Buras et al. (2005b) for ini- 
tial main sequence stars of 8.8 (Nomoto 1984, 1987) 
and 11.2 Mq (Woosley, Heger, & Weaver 2002). Echoing 
the properties of AIC progenitors, the low envelope mass 
and fast declining density (and, therefore, accretion rate) 
are key beneficial components for the success of neutrino- 
driven explosions, but the same properties are also why 
the explosion is necessarily underenergetic. 

The various curves also show a few dips and bumps. 
The first bump, most pronounced in the 1.92-Mq model, 
is associated with an early outflow that eventually fell 
back to smaller radii, while subsequent bumps are caused 
by episodic mass loading of the neutrino-driven wind, 
which sets in ~200ms after core bounce and drives a 
8x1O"3M0 s"^ (5x1O"3M0 s"^) mass loss rate in the 
1.46-M0 (1.92-M0 ; see also Fig. [THl top panel) model. 
The higher-mass flux in the 1.46-M0 model results from 
the higher neutrino luminosity, higher mean neutrino 
energies, and bigger opening angle of escape for the 
neutrino-driven wind. Interestingly, this mass flux is 
strongly angle-dependent, varying by a factor of a few be- 
tween the pole and the angle that grazes the pole-facing 
side of the disk. As described in 21 the dynamical effects 
of the neutrino-driven wind are to entrain the material 
lying along this interface, tearing the disk via Kelvin- 
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Fig. 16. — Color maps of the total neutrino flux soon after bounce (top row) and at the last simulated time (bottom row), for the 
1.46-M0 model (left column) and 1.92-Mq model (right column), with velocity vectors and density contours (for every decade starting 
from 10^* gcm~^). Note the growing flux anisotropy with time, more manifest for the 1.92-M0 model, with regions of high flux coinciding 
with regions of low density, following pri marily the angle-dependent surface area of the neutrinosphere and, secondarily, the latitudinal 
dependence of the temperature (Fig. llUi near such neutrinosphercs. Note the gravity darkening effect d la von Zeipel. From top to 
bottom, left to right, the vector of maximum length corresponds to a velocity magnitude of 52500 km s^^ (infall), 53000 km s^^ (infall), 
18500 km s~l (outflow), and 15500 km s^^ (outflow). 



Helmholtz shear instabilities and mass-loading the wind 
along the corresponding latitudes. Further in, neither 
this wind nor the neutrinos have an appreciable dynam- 
ical impact in driving the disk material outwards, a fea- 
ture only excacerbated by the reduced neutrino flux at 
low latitudes. In Fig. we show in the bottom panel 
the latitudinal variation of the asymptotic velocity (solid 



line) and the density (dotted line). Both show an over- 
all decrease towards lower latitudes by a factor of three. 
The dip in density and higher values of the velocity along 
the pole to around 70° latitude are possibly due to wind 
mass loading, in combination with centrifugal support at 
the neutrinosphere for off-polar latitudes. In Fig. E| we 
show at early times (top row) and at the last simulated 
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Fig. 17. — Color map of the density, in the inner ~70km, overplotted with white vectors representing the energy-integrated electron- 
neutrino flux at the last time simulated for the 1.92-Mq model, presented for comparison with the bottom-right panel of Fig. 1161 Note 
how the flux vectors tend to point in a direction perpendicular to the local isochores, modulo variations, e.g., due to the non-uniform 
temperature distribution, but in general conformity with the von Zeipel theorem which stipulates that the radiative flux should be aligned 
with the local effective gravity. The vector of maximum length corresponds to a flux mag nitudc of 3x10^'^ erg cm ^ s ^ . For a more 
quantitative assessment of the electron-neutrino flux magnitude, see Fig. inland Fig.lTOl 



time (bottom row) for the 1.46-M0 model (left column) 
and the 1.92-M0 model (right column) color maps of the 
total neutrino flux in the radial direction, with isoden- 
sity contours overplotted as white curves, and velocity 
vectors as black arrows. First, due to the history of the 
collapse, the neutron star is relatively devoid of overlying 
material in the polar direction, while for the higher-mass 
progenitor a massive (^O.BM©), dense (10^~^° gcm~'^ 
), near-Keplerian disk obstructs the neutron star at lat- 
itudes <±40°. Given this configuration, conditioned es- 
sentially by the mass distribution of the progenitor white 
dwarf, the dynamical effect of a spherically-symmetric 
neutrino flux would be enhanced along the "excavated" 
polar direction. Indeed, we see a strong neutrino-driven 
wind in the polar direction that does not exist in di- 
rections within ^ ±40° of the equator. However, even 



in the absence of this anisotropic matter distribution, 
Fig-Elreveals the strong latitudinal variation of the neu- 
trino flux at a given Eulerian radius, a variation that 
is established independently of the configuration of the 
circum-neutron star disk material. What controls the 
flux geometry is the combination of two effects. First, the 
exceptional elongation of the neutrinospheres along the 
equatorial direction leads to a decoupling radius (surface) 
about 10 (100) times bigger for a polar observer than 
for an equatorial observer. The angle-dependent decou- 
pling radius of neutrinos mitigates this result (Walder 
et al. 2005), but, as shown in Fig. ^1 the latitudinal 
variation along different directions persists in the total 
neutrino flux. Similarly, Fig. 1171 shows the anisotropy 
of the Vf, neutrino flux, rendered by the corresponding 
flux vectors. Notice how the base of the flux vectors in 



22 









U-1 


.46 


Ma 






1.5 


- 








- 




1-0 
























M=1.92U© 


M [10- 
















0.5 














0.0 













0.0 0.1 0.2 0.3 0.4 0.5 

Electron Fraction 

Fig. 18. — Distribution of the total outflow mass (cells with 
outward radial velocity of at least 10000 kms~^) as a function of 
electron fraction, taken at the last time for both models (solid line: 
1.92-M0 model; broken line: 1.46-M0 model). The total mass in 
the outflow, making allowance for any mass lost through the outer 
grid radius, is 0.003 Mq for the 1.92-M0 model and 0.004 M© for 
the 1.46-Mq model. 



the high-density central regions is perpendicular to the 
local isodensity (or, equivalently, equipotential) contour. 
Second, as shown in Fig. 1101 the temperature and its 
radial-gradient along a given isodensity contour are both 
significantly lower along the equatorial direction, leading 
to reduced diffusive fluxes. These properties are reminis- 
cent of the effect of gravity darkening (von Zeipel 1924) 
in fast rotating (non-compact) stars and the associated 
scaling of the radiative flux with the local effective grav- 
ity (see Owocki et al. 1996), although this may be the 
first time it is reported in the context of a protoneu- 
tron star (but see Walder et al. 2005). With such a 
polar-enhanced wind, the angular momentum loss rate 
is reduced, with consequences for the spin evolution of 
the PNS. 

9. EJECTA COMPOSITION 

In Fig. ^1 we show for the two baseline models the 
electron fraction (Yc) distribution of the material in the 
ejecta (material outside the neutron star moving out- 
wards with a radial velocity greater than 10000 kms^^), 
accounting as well for the mass loss through the outer 
grid radius. Such cumulative outflow amounts to 
4x 10-3 Mq for the 1.46-Mo model and SxlO^^^^foj. 
the 1.92-M0 model. We obtain double-peak profiles, the 
first blast propelling symmetric material {Yc = 0.5), sub- 
sequently followed after 200 ms by progressively neutron- 
rich material, i.e., Y^, = 0.25 — 0.35, in the neutrino-driven 
wind. Note that for these runs we enforced an upper limit 
of 0.5 to the computed Y^ values. Fryer et al. (1999) 
obtained an ejecta mass in the vicinity of 0.2 Mq , two 
orders of magnitude larger than our values. Because our 
ejecta masses are much smaller, we find that the mass 
loss rates and the kinetic energies associated with the 
neutrino-driven wind are relatively more important for 
the global energetics of the AIC of white dwarfs. 

At late times, the asymptotic electron fraction Y^ of 
the neutrino-driven wind varies with latitude (despite the 
smooth variation of other quantities at correspondingly 
larger distances). Material ejected within 20° of the pole 
has an electron fraction of ~0.5, while towards the equa- 
tor, this electron fraction decreases to 0.3, rising again to 



near 0.5 values in regions belonging to the disk (see bot- 
tom row panels in Fig.^. The Y^ values seen in our sim- 
ulations are in fact already set when wind material leaves 
the vicinity of the neutrinosphere, whose properties de- 
pend on the particle trajectory under scrutiny. We know 
from previous studies (Qian & Woosley 1996; Wheeler 
et al. 1998; Thompson et al. 2001; Pruet et al. 2005; 
Frolich et al. 2005) that the asymptotic electron frac- 
tion of the ejecta is controled by competing factors. The 
electron and anti-electron neutrino luminosities, modu- 
lated by the hardness of their respective energy distribu- 
tions, influence the electron flavor production rates via 
the reactions l'el^—^pe~ and UeP^e'^n and, thereby, the 
neutron-richness of the ejecta. The expansion timescale 
sets the duration over which interactions between neu- 
trinos and nucleons can take place. 

The starting value of the electron fraction, i.e., at the 
base of the outflow, is altered by the above factors and 
differs from the asymptotic value seen. In Fig. ^1 we 
show a color map of the electron fraction in the inner 
200 km, highlighting the butterfly shape of the delep- 
tonized region in cross section, in stark contrast with the 
corresponding near-spherical shape seen in core-collapse 
simulations of both rotating and non-rotating progen- 
itors (Keil et al. 1996; Walder et al. 2005; Dessart 
et al. 2005). Deleptonization obtains preferentially in 
the vicinity of the dumbbell-shaped neutrinosphere, and 
stretches outwards for off-polar latitudes. Along the 
equator, deleptonization ceases at smaller radii due to 
the lower effective temperatures (tied to the neutrino 
fluxes). Temperature and neutrino flux are in fact inter- 
twined, since energy deposition by neutrinos may raise 
the temperature locally in the so-called gain region. This 
is also vividly represented in Fig. El (right) by the net 
gain associated with electron-neutrino energy deposition 
in this inner region, which also shows the same butter- 
fly shape. Electron neutrinos emerge from the neutron 
star, and, due to the dumbbell neutrinosphere morphol- 
ogy, at much smaller radii along the poles than for off- 
polar latitudes. The decreasing neutrino flux (dilution) 
reduces this energy deposition beyond ^50 km, and even 
at smaller radii along the equator due to the additonal 
flux reduction there (Fig. I16|l. 

The asymptotic value of the material electron fraction 
is determined in the vicinity of the neutrinosphere and, 
therefore, is directly influenced by this configuration of 
the inner Yc distribution. Along the pole, the wind car- 
ries initially low Y^ material that absorbs electron neutri- 
nos, whose associated luminosity is one magnitude higher 
than that of the anti-electron neutrinos (see Fig. ^ , rais- 
ing the electron fraction to the ceiling value of 0.5 ar- 
tificially adopted in these calculations. Away from the 
pole, the neutrinosphere is located further out, and de- 
spite similar neutrinosphere Yg values, the larger distance 
from the neutron star implies a reduced electron-neutrino 
luminosity and a reduced absorption of neutrinos, lead- 
ing to asymptotic values of the electron fraction of only 
~0.25, not far from the values at the corresponding neu- 
trinosphere. To summarize, the progressive decrease of 
the electron fraction (and of the entropy) away from the 
pole is a result of the reduced electron-neutrino lumi- 
nosities in the vicinity of the latitudinal-dependent neu- 
trinosphere radius (and the associated reduced heating 
and electron-capture rates). 
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Fig. 19. — Left: Color map of the electron fraction in the inner regions of the 1.92-Mq model at the last computed time, with isodensity 
contours (white lines; every decade, starting at lO^^gcm""^) and velocity vectors (black arrows) overplotted, showing in particular the 
locations of the deleptonized regions and the neutrinospheres. Notice the Yc change in the material as it moves away from the neutrinosphere 
(roughly the 10^""^^ gcm^'^ contour). It rises steeply from 0.05 to 0.5 in just a few tens of kilometers along the pole, while this rise is more 
gradual at ~45° latitude, saturating at ^0.25. Right: Same as left, but for the net energy gain due to absorption (positive) and emission 
(negative) of electron neutrinos (the colorbar range used saturates the gain values in the inner regions). Notice how much stronger (and 
positive) the gain is along the poles, causing the Yc of the material to rise steeply to 0.5, while its more modest value away from the pole 
alters the Yc of the corresponding material only up to ~0.25. Note that the anti-electron neutrino flux is one order of magnitude lower 
than that of the electron neutrinos, and thus, too weak to cause the reverse process of a decrease of the Yg of the material. The vector of 
maximum length corresponds to a velocity magnitude of 9130 kms~^. 



10. GRAVITATIONAL WAVE SIGNATURE 

We estimate the gravitational wave emission from as- 
pherical mass motions in our models via the Newtonian 
quadrupole formalism as described in Monchmeier et al. 
(1991). In addition, we compute the gravitational wave 
strain from anisotropic neutrino emission employing the 
formalism introduced by Epstein (1978) and developed 
by Burrows & Hayes (1996) and Miiller & Janka (1997). 

Due to rapid rotation and the resulting oblateness 
of the core, one would expect that rotating AIC mod- 
els would have significant gravitational wave signatures. 
However, though the contribution to the metric strain 
in the equatorial plane, of the aspherical and dy- 
namical matter distributions is not small, we find that 
that of the aspherical neutrino field is larger in mag- 
nitude, though at much smaller frequencies. While we 
calculate that /i+(max) for the matter in the 1.46-Mo 
model is ^--^5.9x10"^^, with a spectrum that peaks at 
~430 Hz, the corresponding /i+(max) due to neutrinos is 
^--^4.6x10"^-^ (derived from the fluxes at 200 km), but at 
frequencies between (~0.1-10 Hz). The total energy ra- 
diated in gravitational waves is ~5.7x 10~^° Mq . The 
corresponding numbers for the faster rotating and more 
massive 1.92-M0 model are ^--^3.6x10"^^ (matter), ~165 
Hz, ~2.0x 10-20 (neutrinos at 300 km), and -7.0x10^* 
M© (? . Note that almost all of the energy is being emit- 
ted by mass motions (99.8% in the 1.46-M0 model and 
98.4% in the 1.92-M0 model), since the power scales 
with the time derivative of the wave strain which is 
small for the waves emitted from the aspherical neutrino 
field. 

We compare the above numbers with those obtained 



by Fryer, Holz & Hughes (2002) for an AIC model of 
Fryer et al. (1999) which was setup with a simple, solid- 
body rotation law and had a final T/|W| of ~0.06 (for 
our 1.46-Mo model: 0.059; see Table 1). They did not 
consider anisotropic neutrino emission. Our more realis- 
tic initial models yield maximum (matter) gravitational 
wave strains that are 1.5 to 2 orders of magnitude smaller 
than those predicted by Fryer, Holz & Hughes (2002). 
The total energy emissions match within a factor of two 
since our models emit at higher frequencies. 

Based on our results, we surmise that gravitational 
waves from axisymmetric AIC events may be detected 
by current LIGO-class detectors if occurring anywhere 
in the Milky Way, but not out to megaparsec distances 
as suggested by Fryer, Holz & Hughes (2002). It is, how- 
ever, likely (§5) that at least the 1.92-JVI0 model will un- 
dergo a dynamical rotational instability leading to non- 
axisymmetric deformation (which can not be captured by 
our 2D approach) and, hence, to copious gravitational 
wave emission over many rotation periods, greatly en- 
hancing detectability. 

11. DISCUSSION AND CONCLUSIONS 

We have presented a radiation/hydrodynamic study 
with the code VULCAN/2D of the collapse and post- 
bounce evolution of massive rotating high-central-density 
white dwarfs, starting from physically consistent 2D rota- 
tional equilibrium configurations (Yoon & Langer 2005). 
The main results of this study are: 

• The AIC of white dwarfs leads to a successful ex- 
plosion with modest energy <10^° erg, thus compa- 
rable to the energies obtained through the collapse 
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of 0/Nc/Mg core of stars with ^8-11 M0 main se- 
quence mass (Kitama et al. 2005; Bm'as ot al. 
2005b). This is, however, underenergetic, by a 
factor of about ten, compared with the inferred 
value for the core collapse of more massive progen- 
itors leading to Type II Plateau supernovae. Al- 
though less and less likely to be the engine behind 
most core-collapse supernova explosions, the neu- 
trino mechanism can successfully power explosions 
of low-mass progenitors and AICs due to the lim- 
ited mantle mass and steeply declining accretion 
rate. 

• Due to high-mass and angular-momentum accre- 
tion, white dwarf progenitors leading to AIC can 
have masses of up to ^2 M© and rotate fast, with 
rotational to gravitational energy ratios of up to 
a few percent prior to collapse. The asphericity 
of such white dwarfs allows the shock generated 
at core bounce to escape along the poles in just 
a few tens of milliseconds, opening a hole in the 
white dwarf along the poles. The blast expands and 
wraps around the progenitor and escapes the grid 
(^5000 km) within a few hundred milliseconds, at 
which time a neutrino-driven wind has grown in 
the pole-excavated region of the white dwarf. Both 
the original blast and the wind show strong lat- 
itudinal variations, partly constrained by the ob- 
structing uncollapscd equatorial disk regions of the 
progenitor, whose centrifugal support prevents it 
from collapsing on a dynamical timescale. Rota- 
tion in such progenitors, thus, affects both directly 
and indirectly the morphology of the explosion. 

• The neutron stars formed have masses on the order 
of 1.4 M0 , with rotation periods close to a millisec- 
ond in the rigidly rotating inner ^^30 km. The fi- 
nal rotational to gravitational energy ratios, for our 
two test cases, cover 0.06 to 0.26, the latter being 
large enough to grow non-axisymmetric instabili- 
ties. At the end of our simulations, the neutron 
stars are oblate and pinched along the poles, with 
polar and equatorial radii in the ratio 1:15 for the 
faster rotating (I.92-M0) model. 

• The morphology of the neutron star leads to a lat- 
itudinal variation of the neutrino flux, the net en- 
ergy gain, and the temperature. In the faster ro- 
tating model, the 'v^" neutrino flux is reduced, 
while the anti-electron neutrino flux is a factor of 
ten lower than that of the electron- neutrino. This 
raises the electron fraction of the ejected mate- 
rial to values close to 0.5 along the poles, but to 
only ~0.25 at lower latitudes, since the correspond- 
ing neutrinosphere is more remote and, thus, the 
electron-neutrino flux is smaller. This introduces 
a latitudinal dependence of the electron fraction of 
the ejected material, but more importantly allows 
neutron-rich material, with entropy on the order of 
20-40 fce/baryon, to escape. Thus, a low-entropy 
r-process might take place under these conditions. 



• The high original angular momentum of the pro- 
genitor follows the mass and is, thus, found mostly 
in the neutron star at the end of the simulation. 
However, rotational energy is also given to the 
ejecta, which uses it to gain (planar) radial kinetic 
energy to escape the potential well, while the rest is 
found in a quasi-Keplerian disk of up to 0.5 Mq in 
the 1.92-M0 model. This disk is an essential com- 
ponent of these AICs; it coUimates the explosion 
and the neutrino- driven wind and also suggests a 
second stage of long-term accretion onto the com- 
pact remnant. 

• The total ejected mass is only of the order of a 
few times 0.001 M0, with only a quarter in the 
form of ^^Ni. The original blast and the short- 
lived neutrino-driven wind will lead to a consider- 
able brightening of the object, but the small ejecta 
mass will quickly become optically-thin, gamma 
rays leaking out rather than depositing their energy 
to power a durable light curve. Therefore, these ex- 
plosions should be underluminous and very short 
lived. Their appearance may also vary consider- 
ably with viewing angle, depending on the mass of 
the progenitor and the presence of a sizable disk in 
the equatorial regions. 

This study has shown that more consistent, rotat- 
ing 2D models alter considerably our understanding 

of accretion-induced collapse, previously obtained un- 
der the simplifying assumptions of spherical symmetry 
and/or zero rotation. Further improvements will come 
by including a consistent temperature structure for the 
progenitor white dwarf and by accounting for the effects 
of magnetic fields. Due to the rapid differential rotation, 
magnetic fields could be amplified considerably and re- 
sult in MHD jets that might alter yet again our overall 
picture of accretion-induced collapse and the energetics 
of the phenomenon. Three-dimensional effects may also 
alter the protoncutron star properties presented here, 
since the faster rotating (1.92-M0 ) model is expected 
to experience non-ax;isymmetric instabilities. 
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